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INTRODUCTION 


The recent thrust in the study of elastic wave scattering has been 
highly motivated by its applications in various fields such as seismic 
explorations, nondestructive testing, material property evaluation and 
dynamic stress concentration. An excellent account of history and 
fundamentals of elastic wave motions is given in [1-2]^ and a comprehensive 
discussion of applications in nondestructive evaluation from a .theoretical 
viewpoint can be found in [3-4] . 

The scattering of a single ellipsoidal scatterer in an infinite 
medium is of fundamental importance and is attacked by the method of 
matched asymtotic expansions [5], the direct volume integral formulation [6], 
the surface integral formulations [7], and also, recently, the polarization 
approach [8] . 

It appears that Mai and Knopoff [6] were first in presenting a direct 
volume integral formulation where they gave the scattered displacements in 
terms of volume integrals involving the displacements and strains inside 
the scatterer. Not knowing these fields, they used the solution when no 
scatterer was present to obtain an approximate solution for a perfect 
sphere. The same approach was later taken by Gubematis [9] for an 
ellipsoidal inhomogeneity. These solutions are appropriate at longwave 
scattering, ka < 1, where k is wave number and ”a" is a typical geometric 
dimension. 

Using the equivalent inclusion method, Eshelby {10,11,12] studied 
the static edastic fields ’’inside" and "outside" an ellipsoidal inclusion 


1 Number in brackets designate References at the end of paper. 




or inhomogeneity embedded in an infinite isotropic elastic medium under 
applied tension. All through this paper, an ’’inhomogeneity’' is referred 
to as a region of different elastic moduli and density compared with its 
surrounding matrix and an "inclusion” is referred to as a region with the 
same elastic moduli and density as its surrounding matrix but include in 
it a distribution of eigenstrains. Eigenstrains are strains that are not 
derived from mechanical loading. As examples we note that thermal strains 
and also the swelling strains due to the presence of moisture are special 
types of eigenstrains. 

The method of equivalent inclusion is a method where the inhomogeneity 
is replaced by an inclusion such that solutions for the two problems are 
exactly the same, the basic concept is sketched in Figure 1. Wheeler 
and Mura [13] first developed but did not apply a complete eigenstrain 
formulation to the dynamic case. The purpose of this study is to extend 
the method of equivalent inclusion to fully develop the equivalence 
conditions and to present a method for complete determination of the 
eigenstrains and/or their derivatives as appropriate. 

The equivalence conditions and the solutions to the scattering of an 
inhomogeneity in terms of eigenstrains and/or their derivatives are first 
developed. Agreement with other approaches is then easily seen. The 
formulation is general and both the inhomogeneity and the host medium 
can be anisotropic. The scattering of an ellipsoidal inhomogeneity in a 
linear elastic isotropic whole space subjected to plane time-harmonic wave 
is studied and the differential and total cross-sections for a uniformly 
distributed eigenstrain are detailed and shown. 

Mr. H.S. Jing and Mr. Y. Paul Hsu assisted in the computer 
programming and graph plotting. Numerical work was conducted by using 
the AMDAHL 470 computing system and plotting was done on the MINC System 
of the Department of Engineering Mechanics, Ohio State University. 
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EQUIVALENCE CONDITIONS 


In this section, we develop the equivalence conditions by requiring 
that the displacement and stress fields obtained in the inclusion problem 
be identical to those for the inhomogeneity problem. Fig 1(a) and 1(b). 

The Inhomogeneity Problem (Problem I) 

Consider the problem of a single inhomogeneity occupying the region 51 
in the infinitely extended region D-51 subjected to applied incident wave 
field > Figure 1. Let the elastic moduli and the mass density be 
denoted by and p' for the inhomogeneity, and by and p for the 

host medium, respectively. 

The governing equations for the displacement field are: 


C., u , + AC., u , = pu. + Ap u. 
jkrs r,sk jkrs r,sk j j 


in D 


in which we used 


a., = C., u +AC., u 
jk jkrs r,s jkrs r,s 


( 1 ) 


( 2 ) 


Ap 




in D-fl 
in 51 


C3] 


AC 


jkrs 


0 

C ' - C 

jkrs jkrs 


in D-51 
in 51 


C4) 


Let the superscripts (i) and (m) denote fields associated with the 
incident wave and the mis-match in mass density and elastic moduli. It is 
clear that 


o 


u . 
J 


(i) ^ Cm) 
u. + u: 

J J 


C5) 


as in the absence of mis-match, disappears and the total field is 
identical to the incident field. 

The boundary conditions are that (1) the displacements and tractions 
at the intersection of the regions and D-J2 must be continuous, and (2) 
the characteristics of out-going wave field and that the stresses die out 
at infinity must be observed. 


The Inclusion Problem (Problem II) 

Consider next an infinite elastic solid of homogeneous moduli 
and density p with distributed eigenstrains , denoted by ^ region $1, 

such that n is identical in shape and size to that in Problem I, and 


e* . = 
ij 



in D-J2 
in 


( 6 ) 


The total strain field is 


e 

rs 


(u + u )/2 
^ r,s s,r*^ 


e® + e* 
rs rs 


where 


a 


jk 


' jkrs 


e 

£ 

rs 


(7) 


C8) 


Using Equations (6,7) in the equations of motion, we easily obtain the 
governing equations for the total displacement field as follows: 


C u . = p u. + C., e* , in D 

jkrs r,sk j jkrs rs,k 


( 9 ) 
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It is clear from Equation (9) that 

u. = (10) 

3 3 3 

where u? is the displacement field due to the presence of and it 

disappears when vanish. The only boundary conditions are those regular 
conditions at infinity and the radiation condition. 

Equivalence Conditions 

For a complete equivalence between Problem I and Problem II, we 
require that the displacement and stress fields in the two problems be 
identical. Hence, for equivalence in stress field, we require, from 
Equation (2) and Equations (7,8), 

^ ^jkrs ^r,s ^ ^^jkrs ^r,s^j ” ^^jkrs ^^^r,s " ^rs^ ^ 

For equivalence in displacement fields we require that the Equations (1) 
and (9) be identical, hence 


Ap u. - AC., u , 
j jkrs r,sk 


'jkrs rs,k 


II 


It is clear that Equations (11,12) are automatically satisfied 
D-S2 by observing the definitions given in Equations (3,4,6). 
to split the RHS of Equations (12) into two parts such that 


C. 


= . 


jkrs rs,k 


^^jkrs ^r,sk 


( 2 ) _ 


‘jkrs rs,k 


= Ap u. 


( 12 ) 

in the region 
It is convenient 

(13a) 

(13b) 
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Employing Equations (11-13) and Equations (5,10) we obtain the equivalence 
conditions as : 


AC., 

jkrs 


r.s'- 


C-i, 

jkrs 


rs 


^^^(r) = -AC. 


jkrs 


u*^^^ (r) 
r,s 


in Q 


(14a) 


Ap (r) + e*^“^(r)=-Ap (r) in 0, (14b) 

These conditions can be used to determine the eigenstrain distribution that 

is necessary for the equivalence of Problems I and II provided that we can 

write u^ ^ in terms of the eigenstrains . One such method is given in [15]. 

It is of interest to note that Equations (14a) are identical in form as the 

equivalence conditions in the static case and that only the components 

of C., are needed for determining Further discussion on the 

jkrs rs,k ® j 

determination of the eigenstrains will follow in the next section. 
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THE SCATTERING OF AN INHOMOGENEITY 


Instead o£ finding the solution to the physical problem stated we 
seek the solution to the equivalent inclusion problem. The governing 
equations are Equations (6, 7, 8, 9). Let the incident wave field be plane 
time-harmonic then the time harmonic displacement and eigenstrain fields 
can be written as, for example, 

u^^'^Cr.t) s uj'™^ (r) exp (-icot) 


C-iwt) 

where w is the frequency of the incident wave field and i*i =-l. Using the 
dynamic version of the Betti-Rayleigh reciprocal theorem and suppressing 
the time dependence we obtain the displacement field as: 


u Cr') = 
m 


u^^^ (f) 

m 


Iff 

V 


'jkrs 




'rs ,k 


(?) dV 


(15) 


or, upon employing Gauss' theorem. 


u (r') = u^^^(r') + /// C g. (r,r’) e* (r) dV 
™ m jkrs ®jm,k 


rs 


( 16 ) 


where g^^(r,r') is the spatial part of the solution to the associated 
Green's function problem for Equation (9). Note that the use of the Green's 
function preserves the characteristics of an outward propagating wave and 
satisfies the boundary condition on stress at infinity. Since there are 
two types of eigenstrains , we write the displacements as 
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C173 




/// Skrs 


- /// 

Q. 


g. , , (r,r*) C., (r) dV 

° 3 m,k’^ '' jkrs rs ^ 


where Equations (6,15,16) are used. It is clear that only the component 
( 2 ) 

of C., e* 1 are needed to determine the displacements. We can therefore 
jkrs rs,k ^ 

view this as a vector quantity, say where 

j 


TT^ 

J 


C., e*^^l 

jkrs rs,k 


(18) 


Substituting Equation (18) in Equation (17) we obtain the solution form 
given as Equation (2.19) in Reference [8], p. 291, for the time-harmonic 
case. If the Equations (14) are substituted in Equation (17), we find the 
solution form Equation (12) in Reference [6], p. 379, or Equation (2.25) 
in Reference [14], p. 2806. 

The stress field that is arisen from the presence of mis-match or 
equivalently by the presence of eigenstrains can be obtained by using Hooke's 
law and Equations (5, 10, 17) as: 


(m) 

pq 



pqmn 






2 ^pqmn ^®jra,k'n' * ^jn,k'm'^ *“jkrs 


( 19 ) 




Quantities of interest such as the differential cross section 
dP(u)/dJ^ defined as: [16,14] 
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( 20 ) 


2 . s *s 

^ ^ 2,. a. . u.> 

dP (gj) _ lim i ii i 

dfi ~ <i°> 


can be obtained in terms of the eigenstrains via Equations (17-20), 


where 1° is the incident power. 


l“ = off 

1 13 3 


<f(t)> denotes time averaging of a function 


<f(t)> - ^ / f (t)dt , 
o 


and n^ are direction cosines for r and n, respectively. 

The differential dfi is the differential element of a solid angle. The 
total cross section is simply 

PW = / ^ (21) 

47T 

In what follows we give examples for the scattering of an isotropic 
smooth inhomogeneity in a linear elastic isotropic infinite medium subjected 
to plane time harmonic incident wave field. 
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LINEAR ELASTIC ISOTROPIC MEDIUM 


For such a medium the spatial part of the Green's function is 

g. (r-f') = . ^ 2 iSR)/R 

4'rrp(o'‘ jm ^ 

+ [(exp i0R)/R - (exp iaRj/R],^^} 

where 


(223 


R = 




and X, u, Vrj. are the Lamp's constants, longitudinal wave speed, transverse 
wave speed, respectively. Before we substitute Equation (22) in Equation (173 > 
we expand the unknown quantities associated with the eigenstrains in form 
of a polynomial as [15,17]: 


TT*(r3 = A. + A, 


J 




^jkJl 


(233 


=‘B.. * x^ * 


(24) 


where A^ , ^ ••• constants. Substituting Equations 

(22,23,243 , with ir^ defined by Equation (183, we obtain 


A. . A.^ * ... 


+ F . . (r3 B. . + F . (r3 B. + . . 

mij ^ ^ ij mijk^ ^ ijk 


(25a3 
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where 


47rpo) f . (r) * 

lUJ 


- 0 <P 6 . * ipf . - <pf . 


2 - 2 

4TTP0J £ ., fr) = - s d), 6 , + U;, . - <t>, 

mjk^ ' mj ^k,mj ^k,mj 


'*’'0“^ = - [Xa^ S.. * 2u ♦,. 6^. 


- 2u . . + 2y (I) , . . ] 


= - [X S. . * 2u S 


mj 


and 


- 2y t|^, . . + 2y *. . . 1 


(|»(r) = /// (exp i6R)/R dV^ 
Q 


= /// (exp i0R)/R dV’ 


*^kJl...s^^^ \ h "' i 3 R)/R dV 


>^Cr) = /// (exp iaR)/R dV 


= /// (exp iaR)/R dV 


C25bj 

(25c) 

(2Sd) 

(25e) 

(26a) 

(26b) 

(26c) 

(26d) 

(26e) 
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(26f) 


= /// ••• (exp iaR)/R dV 


The ({)- and i|;>integrals given in Equations (26) are the volume integrals 

associated with the inhomogeneous Helmholtz equation. They can be carried 

out for an ellipsoidal region by expanding (exp ikR)/R in Taylor series 

expansions with respect to r', for r > r’ and with respect to r for r < r'. 

Here k can be either a or 6. Details are given in Reference [18]. This 

type of expansion for the integrand is particularly useful in determining 

the coefficients of a ’'polynomial” distribution of tt* and A*.. 

J ij 

To determine the coefficients A., A., , B., , B,, , we 

3 JK )K j Kx. 

substitute Equations (18,23,24,25) in Equations (14) and note that we are 

dealing with time-harmonic displacements. Since the tt* and are given 

in terms of polynomials we expand the u?”*^(r), u^"*^ (r) and u^^^(r) in Taylor 

J ^ j s r , s 

series expansions with respect to the coordinate origin , by matching the 
coefficients of terms to the same power of x? x x. x .... we obtain a 

^ 1 1 1 J * 

set of infinite number of algebraic equations for A., A. ,, ... B.,, B.,, .. 

J 3 ^ J k J kil 

The scattering cross sections are given in terms of these coefficients by 
way of Equations (17,19,20). 
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EXAMPLES: 


Uniformly Distributed Egens trains 


Let the region be an ellipsoidal region of 2a^, 2a2 and 2a^ along 

the X, y, and z-axis, respectively. Figure 2. We first expand the integrals 

and their derivatives in Taylor series for r < r’ and obtain the Taylor 

series for (r) from Equations (25). Substituting this series for uj™^ (r) 

and its derivatives in Equations (14) with Equations (23,24) and the Taylor 

series for u^^^ (r) , we obtain the governing algebraic equations for A., A.,, 

J 1 J 

..., ... by comparing the order in the power series i.e. x?, x^, 

x^x^. , etc. To save space these equations are not shown here. Once these 
coefficients are determined, we go back to Equations (25) and find the 
scattered displacement field which is u?(r) when r -► «. By using formulas 
given in [18] the scattered displacement field is given in terms of a triple 


sum. 

In what follows we consider the case for a plane time harmonic wave 
propagating in the +z -direct ion, i.e. 


Uj^^(r,t) = u^ exp i(az - ojt) 


(27) 


For a given (ka) enough terms must be taken in determining the coefficients 
A's and B's in Equations (23,24) from 

[f] {A} + [F] {B} = {H} 

mXm mXl mXm mXl mXl 

[d] {A} + [D] {B> - {E} 

mXm mXl mXm mXl mXl 

where [f ] , [F] are defined by the Taylor expansion of the f- and F- functions 


13 



at r = 0, and [d] , [D] are the average of the synunetric part of [f] , [F] . 


The RHS are obtained from the Taylor series expansion for u. and e;. 


(i) 


(i) 


ij 


respectively. 


Using Equations (12,17) in [18] and Equation (25), we find that, for 
the lowest order of and rr* or ^rs^^ " 

constant terms in Equations (23,24), at a distance far away froir the 
ellipsoid: 




r -► 00 


« -1 <» n n-S, 

(4TTP0) ) III 
n=0 1=0 k=0 


(- 1 ) 4TTa.^a^a^ ^ 1 ^ 2^3 ^ 1 ^ 2^3 

(2n+3) (2n+l)n: (1/2) (k/2) ! (n-iL-k)/2i 


J exp iar 
r 


r-a^(ia)” Z Z.k. - Xia^(ia)^Jl 5, .B, . 

L Ja 1 1 ^ m kj kj 


2uia^ (iot)^ Z Z.Z. B. . "I 
^ m k j kj J 


+ [S^(Jt Z.~6 .) A. - 2ui6^ Z,6 .B, . 

r ^ m j ray j k mj k] 


+ 2yi8^ Z Z.Z. B ] (i3)^] exp (-iut) 

rakjicj-' ^ 


(28) 


where the repeated subscripts must be summed from 1 to 3. Note that B^^ 

is non-dimensional and is homogeneous in (au ) and A. is of dimension 

o j 

2 

(Ap • 0 ) • u^) . The first term in the expression, i.e. n = 0, gives the 

following: 
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(29a) 


(r) 
m 

(aa 


u ar in'- Sr in'- 

1 o 


where Cr,9,(f) are spherical coordinates and 


G°(0,« = - (a^aj/sa^a^) [i^i .AJ (ap/p) 


* (l-a"/s"3 B*. * 2 BJ.] 


C29b) 


H°C0,<i)) = [(B/a) At CAp/p) 


- 2 »/»)'■ \ » 2(B/a)“ 


(29c) 


in which 


At 

3 


A^/(AP0J u^) 




(29d) 

(29e) 


By using Equations (20,27,28) we obtain the differential scattering 
cross section, after manipulation, as 


=a^(9,<t.) + (a/6) a^(9,<t>) 


(30) 


where 


o^'ce,*) = Coa^)® S“ce,4>) • C^il 
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3 ^( 6 ,*) = Caa^)® • H°Ce,((.) H°C9,4.) CD^ • B„,) 


“ n n-1 C-i)"C“a.)"Ca,/a,)®'Ca /a 

C = I I I ^ 2 13 1 12 3 


n n^o 1^0 kio C2n+3) (2n+l) n! a/2)! (k/2) ! (k/2! ! Cn-Z-k)/2! 


« n n-Jt C-i)"®a )“Ca,/a 

0- I I I 


n=0 5.=0 k=0 


C2n+33 C2n+l)nl (4/2] (k/2) ! (n-4-k3/2; 


The super bar here denotes complex conjugate. The total differential cross- 
section can be easily obtained as 


P‘(»)/aj^ = Co a^)‘‘ / G°i:e,4.) G^Ce,^ dn 

I] 

= (a/S)2(a dfi 
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Discussion 


An eigenstrain approach to the scattering of a single ellipsoidal 
inhomogeneity is studied in detail. A complete formulation of the 
equivalence between the inhomogeneity problem and the inclusion problem 
is given while the case Ap=0 was given in NASA contractor report #3445 [15] . 
It is shown that this approach is identical to other approaches such as 
the direct volume integral formulation. Refs. [6,8,9]. 

The eigenstrains are expanded as a geometric series and the coefficients 
are determined by the equivalence conditions derived here and also Ref. [19] . 
The scattered displacements and stresses are given in an analytic series 
form in terms of these coefficients. The scattering of an ellipsoidal 
inhomogeneity in an isotropic elastic mediian under plane time-harmonic 
incident wave is worked out as an example. The eigenstrains are asstuned 
to be uniformly distributed. 

The advantages in using the approach taken is at least three-fold: 

(i) the radiation condition for out-going waves and the continuity con- 
ditions at the interface between the inhomogeneity and the matrix are 
automatically satisfied, (ii) the solution for special inhomogeneity 
geometric shape, such as sphere^ cylinder, plate, disk, prolate and oblate 
spheroids, can be easily obtained by setting the appropriate ratios between 
a.^, 3.^,3.^ in evaluating the volume integrals [18], (iii) with the assistance 
of asymptotic expansion method the solution can be easily modified to 
obtain solutions appropriate for different ranges of wavelengths that are 
of interest. 
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Numerical results and graphical plots for an inhoraogeneity of 
spheroidal geometry, Fig. 2, are given in Figs. 3-9, for the case of 
uniformly distributed eigenstrains . A comparison of the results obtained 
with "exact" solution can only be done for simple cases. Such is done 
for the cases of a planar or a spherical geometry. The comparison is 
given in a subsequent report. Finally, it should be mentioned that a 
generalization of this work for investigating the attenuation and velocity 
factors is being planned. 
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Fig. 1 (a) The inhomogeneity problem, 

(b) The inclusion problem. 




Fig. 2 An ellipsoidal inhomogeneity under incident wave 
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A2/A 1=2/3, A3/A 1=2/3 
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ALPHA A1 

Fig. 4 G(9,(t)) VS. aa^ for different ratio of a .a a,/a^: 

titanium in tungston, (j) =0, 9 = 0, G*G = 
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Fig. 3 


.4 .8 1.2 1.6 2 

ALPHA A1 


G(0,(^) vs. for different ratios of 

a^/a^: tungston in titanium, (|) = 0,9 * 0, GfG = 
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Fig. 8 Polar plots for the scattering amplitudes for the 
case aax = 2.0, = 0, a2/a^ = 3.^/ a^ = 2/3. 
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Fig. 9 Polar plots for the scattering amplitudes for the 
case aa^ = 2.0, 4 > = 0, ~ ^3^^1 ” 
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